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Abstract. Two models based on the hydrostatic primitive equa- 
CN| ■ tions are proposed. The first model is the primitive equations 

with partial viscosity only, and is oriented towards large-scale wave 
p j structures in the ocean and atmosphere. The second model is the 

' viscous primitive equations with spectral eddy viscosity, and is 

oriented towards turbulent geophysical flows. For both models, 
' the existence and uniqueness of global strong solutions is estab- 

d . lished. For the second model, the convergence of the solutions to 

the solutions of the classical primitive equations as eddy viscosity 
parameters tend to zero is also established. 

>* ' 1. Introduction 

t> 

We study two models for geophysical flows based on the hydrostatic 
^ . primitive equations; both are designed to faithfully simulate certain 

phenomena in the geophysical flows but they are motivated by differ- 
ent physical considerations. A distinctive characteristic of the flows 
under consideration is that the vertical scale (~ 10km) is much smaller 
than the horizontal scale (~ 6000km). Thanks to this disparity, a 
hydrostatic approximation is possible, and gives rise of the primitive- 
^ ; equations. 

■ On the mathematical side, the theory for the primitive equations is 

fairly complete; see, e.g., the pioneering work in [22, 23], and the survey 
article [28]. In particular, in contrast with the Navier-Stokes equations 
[9, 33], the primitive equations have been shown to have unique global 
strong solutions [7, 20, 21, 18]. 

The first model we study aims to faithfully simulate large-scale co- 
herent structures including wave phenomena in the ocean and atmo- 
sphere. For relevant discussions of this topic, see, e.g., [12, 24, 27]. 
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For the phenomena that we are interested in, both the ocean and at- 
mosphere are close to being inviscid. Therefore, the inviscid primitive 
equation is the preferred model. However, it is very costly to simu- 
late the inviscid primitive equations directly due to the small scales 
embodied within the model. What we propose here is a new model 
with eddy viscosity added to the small scale (high frequency) part only 
while keeping the large scale (low frequency) features intact (at least 
directly). Intuitively, this type of model would reduce the complex- 
ity due to the damping on small scales whereas keeping the desirable 
large scale structures. Such a naive approach may not work all the 
time due to the cascade of energy induced by the nonlinear advection 
term. However, we can easily find situations where such cascade is 
small or negligible. Indeed, it is easy to find exact large-scale solutions 
to this primitive equation with partial viscosity; see below. Of course, 
the existence of such examples do not fully justify the model, and ex- 
tensive numerical experiments are called for which is our future plan. 
The idea of partial damping on the high frequency components of the 
system is not alien to the geophysical community (see, e.g., [11, 25]) 
or the mathematical community (see e.g. [26, 4, 5, 10]), although the 
application to the primitive equations is new here. Our goal in this 
paper is to demonstrate the global well posedness of this model with 
partial viscosity. 

The second model is oriented at turbulence modeling for geophys- 
ical flows. The simulation of three-dimensional turbulent flows is a 
formidable task due to the need to resolve the small scale fluctuations 
or eddies that have subtle effects on the large-scale dynamics of the 
flow. To make this problem computationally tractable, these effects 
must be modeled. In one approach, the velocity field is averaged over a 
small radius to derive equations in terms of the averaged velocity. For 
nonlinear equations, there arises the problem of closure because the 
product operation is not closed under the averaging process. To obtain 
a closed system, the average of the nonlinear terms in the equations 
must be approximated and expressed solely in terms of averaged quan- 
tities. The way in which this is done gives rise to a variety of models. 
The approach we consider, called the eddy-viscosity method, treats the 
Reynolds stress as a viscous effect caused by the transport and dissi- 
pation of energy due to the small-scale eddies. For this reason, this 
additional viscosity is called the eddy viscosity or turbulent viscosity. 
The turbulence model of Smagorinsky [31] belongs to this type. For an 
overall survey on issues related to these models, see [3]. 

Unfortunately, a straightforward application of the approach de- 
scribed above leads to the over smearing of the large-scale structures 
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in the fluid. To remedy this unwanted effect, it has been proposed 
that the eddy viscosity be added only to the subgrid scales. In this 
way, one hopes to prevent the large-scale structure from being smeared 
away. Here, we examine a particular class of models of this type called 
spectral-viscosity or spectral-vanishing-viscosity models, in which the 
scales are defined in terms of Fourier modes. The subgrid viscosity 
is simply realized as an addition of the artificial viscosity only to the 
high-frequency modes. The most intuitive way of doing this to insert a 
high-pass filter to the standard eddy viscosity. This approach was con- 
sidered in [13] for hyperviscosity on the Navier-Stokes equations and 
in [15] for nonlinear as well as hyper-viscosity on the Navier Stokes 
Equations. In both works, the well-posedness of the resulting spectral 
viscosity is proven, and the consistency of these model with the original 
Navier-Stokes equations is discussed. 

We employ the idea of spectral viscosity to build and analyze a tur- 
bulence model for the geophysical flows in the ocean and atmosphere. 
As mentioned above, the primitive equations, even without any eddy 
viscosity, have been shown to have unique global solutions, provided 
the initial and boundary data are sufficiently smooth. For our model, 
we prove its global well posedness, which should not come as a surprise. 
In addition, we will show the convergence of the solutions of the model 
to the solutions of the primitive equations without any eddy viscosity, 
as the eddy viscosity parameters tend to zero. This is not possible 
yet for the Navier-Stokes equations because there convergence is only 
shown to be in a weak sense. 

We should point out that this technique, usually under the name 
of spectral viscosity or spectral-vanishing viscosity, is known in terms 
of turbulence modelling (see [19, 13, 2, 15, 32], and below), and to 
applications in geophysical fluid dynamics [11]. However, to the best 
of our knowledge, the well-posedness result for the three-dimensional 
nonlinear primitive equations with partial viscosity is new. 

The paper is organized as follows. In Section 2, we introduce and 
prove the well posedness of a model with only partial viscosity. In 
Section 3, we introduce the linear spectral eddy- viscosity model. In 
Section 3.2, we prove the existence and uniqueness of strong solutions 
to the model. In Section 3.3 we study the convergence of the solutions 
of the model to the solutions of the original primitive equations. 

2. A MODEL WITH PARTIAL VISCOSITY 

In this section we study a model with partial high-frequency viscosity 
only. The lower modes are not damped directly. This feature renders 
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the model suitable for large scale coherent structures in the ocean and 
atmosphere, because non-physical large scale damping could change 
the large scale coherent structures over time. 
The model reads 



— + (u • V)u + w— + /k x u H Vp - 

f,A(/ - P„,„)u - - P„,„)u = F, (2.1) 

| - (2.2) 

V-u+ — = 0, (2.3) 



In the above, u = (u,v) is the horizontal velocity, w is the vertical 
velocity, \i and v are the horizontal and vertical kinetic viscosities, 
respectively. We use V and A to denote the 2D horizontal gradient and 
Laplacian operators, respectively. The operator [1 — Pm,n) represents 
the high-pass filter and will be defined later on. 

We consider a rectangular domain Q = M x (— H, 0), with M = 
(0, L x ) x (0, L y ). We consider periodic boundary conditions in both the 
x and y directions, and free-slip, non-penetration boundary conditions 
in the vertical directions. More precisely, 



u(x,y,z,t) = u(x + L X7 y t z,t), (2.4) 

u(x,y,z,t)=u(x,y + L v ,z,t), (2.5) 

du . du . . 

7T U=o = 7T \z=-h = 0, (2.6) 
on on 

w\ z=0 = w\ z= - H = 0, (2.7) 

p\z=o = Po- (2.8) 



Under the settings just described, we can define the high-pass filter 
in terms of Fourier frequencies. Specifically, for each function u G 
(L 2 (fi)) 2 we let 



G m ,ne imx 'cosnz (2.9) 
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where 



m = (mi, 1712) G Z 2 , 
|m| sup = maxflmil, |m 2 | 



JUL 1. 



x' = 2tt 

. TTZ 

u m ,n = / u(x, y, z)e~ imv: cos nz' dx dy dz. 
Jn 

Remark 2.1. The issue of suitable physical boundary conditions for 
the inviscid primitive equations is an unresolved one. Partial results 
concerning the linearized primitive equations are available in a series 
of papers [29], [8], and [30]. There, an infinite set of nonlocal bound- 
ary conditions were proposed, which guaranteed the well-posedness of 
the linearized system. Here we avoid this issue and use the periodic 
boundary conditions on the lateral boundaries. 

Remark 2.2. As we have touched upon in Introduction, we can eas- 
ily identify some large-scale motions which are exact solutions of the 
inviscid primitive equations. For example, let 

, . / 47raA / 2ny 
ip — sin — — cos 



Li J \L 

and 

u = 

With the surface pressure po given by 

2tt 2 /4ttxV 8tt 2 /2ny 



(u, 0) is a set of exact solutions of the inviscid primitive equations, 
that is, the system (2.1)-(2.8) without the viscosities, or with partial 
viscosity with M > 4. 

2.1. The barotropic and baroclinic modes. As usual, p and w can 
be expressed in terms of u. Specifically, integrating (2.2) from z to 0, 
and using the boundary condition (2.8) we obtain 

p(x, y, z, t) = p (x, y, t) - p gz, (2.10) 

Integrating (2.3) from z to 0, and using the boundary condition (2.7), 
we obtain 

V-ud£ = V- / ud£. (2.11) 
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Setting z = —H in (2.11), and using (2.7) again, we find 

»o 

udz = 0. (2.12) 

H 

We substitute (2.10) and (2.11) into (2.1), and obtain a single closed 
equation for u: 

— + u ■ V u + V ■ / ude — + /k x u + -Vp - 

dt J z dz p 

yA(I - P m ,n)u - v-q^{I - Pm,n)u = F. (2.13) 

It turns out essential to rewrite this equation and work with the fol- 
lowing form: 



aF + (u ' v)u+ ( v '/ ud «j& +/kxu + ^ Vp »- 

d 2 d 2 
[iAu - v-q^ u = F - /iAP MjAr u - v—P m ,nU- (2.14) 

We note here that Pm,n~u contains only finite number modes of u. The 
prognostic variable u satisfies the equation (2.14), the constraint (2.12), 
and the boundary conditions (2.4)-(2.6). To complete the system, we 
also require u to satisfy the following initial condition: 

u(x, y, z, 0) = u (x, y, z). (2.15) 

The diagnostic variables w and p are given by (2.11) and (2.10) respec- 
tively. 

We now specify the barotropic and baroclinic modes of equation 
(2.14). We let 

1 f° 

u(x,y,t) = —J u(x,y,z,t)dz, (2.16) 

which denotes the barotropic mode of the primitive variables. We also 
denote by 

u'(x, y, z, t) = u(x, y, z, t) - u(x, y, t) (2.17) 
the baroclinic mode of the primitive variables. It is easy to see that 

= 0. (2.18) 
From (2.4) and (2.5) we derive the boundary conditions for u: 

u(x,y,t) = u(x + L x ,y,t), (2.19) 

u(x,y,t) = u(x,y + L y ,t). (2.20) 
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By (2.12) we see that u satisfies the following constraint: 

V-u = 0. (2.21) 

Then u' also satisfy the periodic boundary conditions on the lateral 
boundary, 

u'(x, y, z, t) = u'(x + L x , y, z, t), (2.22) 
u'(x, y, z, t) = u'(x, y + L y , z, t). (2.23) 

It inherits the boundary conditions for u on the top and bottom, 



du' 

dz 



z=0 



du' 

dz 



0. 



(2.24) 



-H 



We now derive the equations that u and u' satisfy by first taking av- 
erage of the equation (2.14): 



du 



du 



1 



d 2 u 



— + (u • V)u + w— + fk x u + — Vp - /iAu - v- 

dt dz p dz 1 



dz' 



/jlAP m ,nU - v—P M)N u. (2.25) 



We notice that 



,d 2 u. 



dz 2 



1 

H 



dz- 



■Pm,nu 





d 2 u 




f 




J-H 


dz 2 




1 




d 2 


~ H 


f 

J-H 


9z2 Pm,n 



0, ( by (2.6)), 



-H 



AP MtN u = AP MN u. 
By using (2.18), we find that 



(u- V)u= (u- V)u+(u'- V)u'. 

Using (2.12), (2.18) and (2.21), we find 



du 



w— = (V • u') ■ u'. 

dz 

Hence the equation for u can be written as 



du 



— + (u • V)u + [(u' ■ V)u' + (V ■ u')u'] 

+ fk x u + — Vp - /iAu 

Po 



F - nAP M , N u. (2.26) 
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The barotropic variable u satisfies the following conditions: 

u(x,y,t) = u(x + L x ,y,t), (2.27) 
u{x,y,t) = u{x,y + L y ,t). (2.28) 
V-u = 0. (2.29) 

Subtracting (2.26) from (2.14) we obtain the equation for the baroclinic 
mode: 



du' 

~dt 



/k x u' - fiAu 



u'-v)u , + (v- J u'eA 

(u' ■ V)u + (u ■ V)u - te • V)u' + (V • u')u' 



du' 

dz 



3 2 u' 

dz 2 

¥'- 



v- 



d 



nAP MiN u' - u—P MtN u'. (2.30) 



In addition u' satisfy 



u'(x, 


y,t) 


u'(x, 


y,t) 


du' 




dz 








f 


u' dz 


J-H 





u'(x + L x ,y,t), 

Ly, t), 



-u{x,y 
du' 

dz 

= 0. 



0. 



-H 



(2.31) 
(2.32) 

(2.33) 
(2.34) 



2.2. Global well— posedness of the model with partial viscosity. 

In order to have uniqueness for the solution of (2.1)-(2.8) we work with 
functions that have zero average over Q. Indeed it can be checked that 
if the initial data and the forcing have zero average over Q, then the 
solutions have zero average over Q at any time. 

We use the convention L 2 (fi), if 1 (Q), etc. to denote function spaces 
that have zero average over Q. We let 

V = {u G C°°(-R 3 ) 2 | u periodic in x with period L x , periodic in y with 

period L y , periodic and even in z ith period 2H}, 

V = V H \ closure of V in (H 1 ^)) 2 ). 

The inner product and norm of H will be denoted as (•, •) and | • |, 
respectively. The space V inherits the inner product and norm of H 1 , 
which will be denoted as ((•, •)) and || • ||, respectively. 
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Since the functions in V has zero spatial averages, we have the fol- 
lowing Poincare inequality for functions in V: 

, du , 



|u| 2 <C7(|Vu| 2 + |^| 2 ) 
oz 



(2.35) 



Therefore (|Vu| 



du 

dz 



2 )2 is equivalent to the usual H 1 norm, and 



can be taken as the norm for V. 

In what follows we abuse the notation by denoting every generic 
constant by C . Such constants may depend on the domain Q and 
the function spaces in the context, but we omit such dependence in the 
notation. But if the constant depends on any other parameters, such as 
M, N etc., we shall use a specific symbol and specify such dependence 
in the notation. 

We shall prove the global existence and uniqueness of strong solutions 
to the system (2.1)-(2.8). 

Theorem 2.1. For a given T > 0, let F £ L 2 (0,T; H), u £ V. Then 
there exists a unique strong solution u £ C([0, T]; V) D L 2 (0, T; H 2 (Q)) 
of the system (2.1) -(2.8) which depends continuously on the initial 
data. 

We shall first obtain some key estimates that will be needed for the 
proof of Theorem 2.1. The proof of the theorem will be furnished at 
the end of this subsection. 



L? estimates 

We multiply (2.14) by u, integrate by parts over Q, using the boundary 
conditions (2.4)-(2.8), we obtain 



[5p u) + ((u ■ V)u, u) + ((V • J° u do|f, u 



+ 



(/k x u, u) H (Vpo, u) + fi\ Vu| 2 + v 

Po 



du 



dz 



(F,u) +ns\VP MtN u\ 



9 P u 

oz 



We notice that 



du 

~dt : 



u 



1 d . 
u 

2 dt l 1 



(Vp , u) n = H(V Po , u) = 0, 
(/k x u, u) = 0. 



(2.36) 

(2.37) 

(2.38) 
(2.39) 
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Let 

6(u, 5, u # ) = ^(u • Vu) + w(u)|^, u#^j , (2.40) 

where w(u) is defined as in (2.11). We can verify that the trilinear 
operator b(u, u, u#) is skew symmetric with respect to the last two 
arguments, that is, 



b(u, u, u*) = — 6(u, u*,u). 
Then it is inferred from (2.41) that 

b(u, u, u) = 0. 
By (2.36), (2.37), (2.38), (2.39) and (2.42) we have 



(2.41) 
(2.42) 



1 d 







|u| 2 + /i|Vu| 2 + v\— u| 2 = (F,u) + 



2 dt 



dz 







v\VP M , N u\ 2 + v\—P MtN u\ 2 . (2.43) 



Since Pm atu contains only finite number of Fourier modes, \VPm n u \ 2 , 
d 

as well as I— Pm n^A 2 can be bounded by |u| 2 . More generally, for 

oz 

functions in H k (Q), there exists a constant Ck(M, N,fl), which is in- 
dependent of the function, such that 



\Pm,nu\ H * < Cfc|ll|L2. 



(2.44) 



Using (2.44) and the Cauchy-Schwarz inequality, we derive from (2.43) 
that 



dt 



u| 2 + dVul 2 + z/|^-u| 2 < CIFI 2 + Cilul 

oz 



(2.45) 



Applying the Gronwall inequality to (2.45) yields 



u(-,t)| 2 + /i / |Vu| 2 ds + ^ 



d_ 

dz 



-u 



ds < J x (t), (2.46) 



where 



Ji(t) = e Clt (\u \ 2 + C |F| 2 ds) . (2.47) 



L 6 estimate on u' 

We take inner product of (2.30) with |u'| 4 u', and integrate by parts 
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over O to obtain 

-—\u'\% + 6(u',u' ,|u / |V)+ ! fkx u'- |u'rVdfi + 
6 dt J n 

J (V ■ V)u + (u ■ V)u' - (u' ■ V)u' + (V ■ uV) |u'| 4 u'dfi - 
J ^Au' + ^j |u / |Vdfi = y F'-lufu'dfi- 

jf (^AP MtN u' + u^P M>N ^ |u'|V dfi (2.48) 

In the above, the trilinear operator b(-, -, •) is defined as in (2.40). We 
can verify by calculations that 



&(u',u', |u'|V) = -2&(u',u', |u'|V), 
which implies that 

6(u', u', |u'| 4 u') = 0. 
Since k x u' is orthogonal to u', we have 

fkxu' ■ |u'| 4 u'dfi = 0. 



(2.49) 



(2.50) 



(2.51) 

Noticing the divergence free condition (2.29) for u and the horizontal 
periodic boundary conditions (2.27)-(2.28) and (2.31)-(2.32) for u and 
u' respectively, we find 

f (u ■ V)u' • |u'|V dQ = 0. (2.52) 
Jn 

For the inner products involving the diffusion terms, we find 

- I ^Au / .|u'|Vdfi = /i / |Vu'| 2 |u'| 4 dfi + /i / |V|u'| 2 | 2 |u'| 2 dfi, 
Jn Jn Jn 

(2.53) 



[ v?^-.\ufu'dn = v / A*\u>\ 4 dn + v [ 
Jn dz 2 J n dz J n 



oz 



|u'| 2 dQ. 
(2.54) 
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For integrals on the right-hand side of (2.48), we use (2.44) (with 
k = 2) to find that 



-//AP/^atu' • |u'| 4 u'df2 

< /ilAP^TvuVlu'Hio 
^Cyu'Lallu' 131 * 



10 



< C|u'|| 2 |u'|J. + \ f |u'| 4 |V U '| 2 dfi + ^ / |u'| 4 
In the above we have used the interpolation inequality 



du' 



dz 



da 



\ L f<miMh, 



(2.55) 



which can be obtained by setting p = 10/3, p\ = 2, p 2 = 6 in (A.l), 
and then using (A. 5). Similarly, 



f ^2 

j^-u—p M>N u'-\ufu'dn 

<C7|u'|i 2 |u'| 4 6 + ^ / |u'| 4 |Vuf &Q + - A [ \u 
4 Jn 4 J n 

Hence we derive from (2.48) that 



/ 1 i 



du' 



dz 



dn. 



1 Vl 6 + 
6d^ |u|i6 + 



J (V • V)u - (u' ■ V)u' + (V ■ u')u') • |u'| 4 u' dfi- 



1 

— V 

2 





|Vu' 


Huf+IVlu'lf 












du' 


2 

|u'| 4 + 


— lu'l 2 

dz 


In ( 


dz 



u'\ 2 I dQ- 



< /'F'-|u / | 4 u / dn + C|u / |i a |u / |i«. (2.56) 
Jn 

For the integrals on the right-hand side of (2.56) that involve nonlinear 
terms, we proceed by integration by parts, using the periodic boundary 
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conditions on u' and u when appropriate, and we find that 

I (u'- V)u- |u'| 4 u'dfi 
Jn 

[ ((v • u')n ■ |u'|V + (u' • v)(|u'| 4 u') • u) an 

Jn 

[ (u- f (V-u')|u'|Vdz + u- / (u'- V)(|u / |V)d^ dxdy 

JM V J-H J-H J 



<C |u| / \Vu\\W\ b dzdxdy. 

JM J-H 

By the Cauchy-Schwarz inequality and the Holder's inequality, we find 



u / |Vu ||u ydzdxdy 

M J-H 



< 



|u| 4 dM 



M 



IVu'Piurdn 





|u'| 6 d^| dM 



M \J-H 



(2.57) 



By the Minkowski integral inequality (A. 6), we have 

»o 



M \J-H 



\ 2 

\u'\ 6 dz\ dM I 



H \JM 



u'| 12 dM dz. (2.58) 



Applying the Ladyzhenskaya inequality (A. 2) to <f) 3 in M 2 , we obtain 



1 12 



Il 12 (m) 



<C\cj>\ 



L 6 (M) 



|4 |V<^| 2 dx dj/ ) + 



M 



1 12 



lL 6 (M)- 



(2.59) 



Using (2.59), we infer from (2.58) that 



M \J-H 



\ 2 

|u'| 6 d^] dM 



< 



< C 



H \ JM 



C / luTdM / |uT|Vu'| 2 dM 



M 



M 
i 



l 

2\ 2 



|u'| 6 dM dz 



H \JM 



<C[ I |u'| 6 dfi 



|u'| 6 dMV (J |u'| 4 |Vuf dM^ + J |u'| 6 dfi 



u'| 4 |Vu'| 2 dfi 



luTdfi. 



14 QINGSHAN CHEN, MAX GUNZBURGER AND XIAOMING WANG 

Therefore we have 



o \ 2 

|u'| 6 d^) dM | 

AI \J-H 



C|u'|| e (/ |u'| 4 |Vu'| 2 dfiJ +|u'||«. (2.60) 
By the Ladyzhenskaya inequality (A. 2) for functions in M 2 , we find that 



'u| 4 dM) <C|u 



M 



(2.61) 



Using (2.60) and (2.61) we infer from (2.57) that 
J (u j |Vu'||u'| 5 d^ dxdy 



< C|u 



3 

,'| a 



L 2 (M) l'^li? 1 (M) I l U lL 6 



uflVufdfi +|u' 



\L<> 



(2.62) 



By Young's inequality, we have 



\ u \h(M)\ u \m(M)\ u '\h [J |u'| 4 |Vu'| 2 dfi 

<C7|u|i 2(A/) |u|^ (M) |u'|| 6 + ^ /" |u'| 4 |Vu'| 2 dfi, (2.63) 

For the other integral that involves nonlinear terms, 

! (u' • V)u' + (V ■ u')u' ■ |u'|Vdfi 
Jn 



M>;.<9 4 (|u'| 4 M ;.)dfi 



Jj J (J < u 'i dz j d A\ u '\ iu 'j) d z) dxd V 
<C [ ( [ \u'\ 2 dz [ |u'| 4 |Vufd^ dxdy. 

J M \J-H J-H J 
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By similar use of the Minkowski inequality (A. 6) and various interpo- 
lation inequalities, we find that 



M \J-H 



\u\ 2 dz I |u'| 4 |VuTdz ) dxdy < 



C|u'|le (n) n |u'| 4 |Vu'| 2 dfi) (|uV (f2) + |Vu'U 2(n) ) . (2.64) 
By Young's inequality, we have 



(|uV (n) + |Vu'| i2( n)) |u'|| 6(n) U |u'| 4 |Vu'| 2 dnj < 

(|u'|| 2( n) + |Vu'|| 2(n) ) |u'|| e(n) + i/* I |u'| 4 |Vu'| 2 da (2.65) 
Putting (2.56)-(2.65) together, we find 

iu'| 6 L6 +/i / (|Vuf |u'| 4 +|V|uf | 2 |u'| 2 ) dfi + 
dt in V / 

u'| 2 ) dfi< (C|F'|| 2 +C 2 |u'|| 2 ) + 



<9u' 



5 I /|2 



9 



C In 



L 2 (M) l U liJ 1 (A/) ^ I^L^M)!"^^/) "I" I" | L 2 (f2) 



U 



.'12 



C|Vu 



'|2 



L 2 {Vt) I l U lL 6 - 



/|6 



(2.66) 



Ignoring the other positive terms on the left hand side of (2.66), and 
dividing both sides by |u'|£ 6 , we have 



dt 



u'|i 6 <C|F'|i 2 + C 2 |u'|£ 2 



C (\u\ 2 L2(u) \u\hm + \"\h { n) + |Vu|| 2(n) ) |u'|f 6 . (2.67) 



Applying the Gronwall inequality to (2.67), and using the L 2 estimate 
result (2.46), we obtain 



W(;t)\ 2 L e < J 6 



(2.68) 



with 



Mt) = eC&W+KiWHW)) (| Uo | ie(n) + c,m +c£ |F|| 2(n) da) . 
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Integrating (2.66) over [0, t], and using the estimate (2.68), we obtain 



H I / I iVuf |u" 4 
'o Jn 



viuflVl 2 ) ands 



JQ, 



du' 



dz 



'|4 



U 



— lu'l 2 

dz 



u'| 2 dQds < J 6 (t), (2.69) 



with 



J 6 (t) = \u \ b L6 + Ji(t)[C I \F(;s)\ 2 L2 ds + C 2 K l tj + 

CJ$(t) (K^+K^t + K^t)) 



Estimate |Vu| L2( - M ) 

We multiply (2.26) by — Au and integrate by parts over M to obtain 



1 d 

2 dt 



Vu|i 2(M) +^|Au|i 2(M) + f /kxu-(-Au)dM 

Jm 



F ■ Au dM — — / Vpo-AudM- f (u-V)u-AudM 

M Po Jm Jm 

- / (u' • V)u' + (V • u')u' ■ AudM + fi / AP^u-Audn. 
Jm Jn 

(2.70) 

We note that, by integration by parts, 



(/k x u) • AudQ 



M 



<\f\ 

|u| L2 |Au| L 2 
4 



<^l/lLluli 2 + ^|Au|i 2 . 



And thanks to (2.21] 



Vp • AudM = 0. 



M 



Following similar steps in the handling of the 2D Navier Stokes equa- 
tions, we obtain 



M 



I 3 

'u- V)u • AudM < C\^l 2 , M) \Vu\ L 2( M) \Au\l 2(M) . 
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Applying the Cauchy-Schwarz and Holder inequalities one have 



/ (u' ■ V)u' + (V ■ u')u' ■ AudM < 

J M 



ClVu'l?,,™ f / |u'| 4 |Vu'| 2 dfi ) 4 |Au 



For this last integral in (2.70), we proceed by the Cauchy-Schwarz 
inequality and (2.44) with k = 2, 



/i I AP M)N u ■ AudQ < /i\AP M>N u\ L 2\Au\ L 2 



< C 2 /i|u| L 2|Au| L 2 

< C 2 n\u\ 2 L 2 + y|Au| 2 2 . 



Then we derive from (2.70) that 

-^|Vu|i 2(M) +//|Au|£ 2(M) < ^|F|| 2(M) + (C 2 ^ + C|/| 2 J|u|i 2 + 

Clul^^lVull.^ + ClVu'H^ + c/ |u'| 4 |Vu'| 2 da (2.71) 

Jn 

Applying the Gronwall inequality to (2.71), and using the previous 
estimates (2.46) and (2.69), we obtain 

|Vu(-,t)|i 2(M) + | jf |Au| 2 2(M) ds < J 2 {t), (2.72) 

with 



Mt) = e CJ ^ (luol 2 ,, + (C 2 pt + C|/| 2 J J ± (t)t + 



C I \Ffa (M) ds + CJi(t) + CJ 6 (t) 



o 



To estimate \d z u\ L 2^ 

We multiply (2.14) by d 2 u/dz 2 and integrate by parts over Q, and we 
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have 



l4:\ d Mh(a) + V\^d z u\ 



2 dt 



L 2 (fl) 



V 



d 2 u 



dz 2 



/ n (,u.V)u + (v./°ud,)g 



L 2 (n) 

du\ d 2 u 



dz 2 



d 2 u 

dz 2 

dn- 



^\VP M:N d z u\l 2{Q) 



dz 2 



Pm,n u 



(2.73) 



L 2 (Q) 



By Holder's inequality, 



d 2 u 

dz 2 



< C\F\ 2 L2 



d 2 u 



dz 2 



L 2 



With regard to the integral in (2.73) that involves the nonlinear con- 
vection terms, we find 

du\ d 2 u 

dz 2 



-/ n (,u.V)u + (v./°„d,)| 

/J( (u .v)u + (v./ua,)g.| d o 



(d z u ■ V)u + (u • V)d z u - (V ■ u 



du 

dz 



V- / ud£- 



d 2 u\ du 



dz 2 J dz 



dtt 



(d z u • V)u - (V ■ u) 



du\ du 



dz J dz 



dfl 



-(V ■ d z u)u ■ d z u — (d z u ■ V)d z u ■ u + 2(u • V)<9 2 u ■ <9 2 udf2. 



n 



Hence 



/ n ( (u .v, u+ (v./ ud ^) 



du\ d 2 u 
) ' ~dl? 



dQ 



< c 



\u\\Vd z u\\d z u\dVL 



< <7|u|£6(n)|9*u|£3(n)|V9«u| L a(n) 

i i 

< C'|uU6 (n) |5 z u|2 2(n) |5 2 u|| 1(n) |V5 z u| i 2 (n) 



< C\u\ L & {n) \d z \Vd z u 



L 2 (Q)\ 



3 

\h(Q) 



< C\u\%\d z u\ 2 L2 + l|u^|| 2 + ^|V<9 2 u|£ 2 
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The last two integrals on the right-hand side of (2.73) can be handled 
by (2.44) with k = 2. Thus we derive from (2.73) that 



dt l 



<^ u Il2(q) + (A Vd z u\ 2 L2(n) + v 



(fi) 



d 2 u 



L 2 (Q) 



<C\F\ 2 L2 



C (Ji(t) + Jl{t)) \d z n\% + C 2 |u|| 2 . (2.74) 
An application of the Gronwall inequality to (2.74) readily yields 

\dzu(-it)\h(n) + V / I V<9 z ii|£ 2(n) ds + v j |— |* a(n) ds < J z (t), 

(2.75) 

with 

J z (t) = e c ^ +J ^y (ju |^ + cj* |F& d + C 2 J^ . 
To estimate |Vu| L 2 

We multiply (2.14) by — Au and integrate by parts over f2, 

5>-li. m + I (<» ■ V)u + (V • f° u d{) I) • (-Au) dfi + 

/ / fcxu-(-Au)dn+- / Vpo-(-Au) df]+/i|Au|^ 2(f7) +z/|V<9 z u|^ 2(c) 
in Po Jn 



(F, -Au) + /i|AP Mi7V u|^2 (n) + v\VP MjN d z VL 



\L 2 (n)- 



(2.76) 



We note that 



and 



/ (/k x u) ■ Audfi = 0, 

/ Vj9 • (-Au) dQ = - [ Vp • / (Au) d^ dM 
Jn Jm J-h 



-H / Vpo-AudM 

<M 



It is easy to see that 



/ o ((».V)» + (v.jf adf)£).<A»)dn< 

C y (ju||Vu| + y |Vu|d^|<9 2 u|^ |Au|da 
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By the Young's inequality and the Sobolev interpolation inequality 
(A. 4) for L 3 functions, we have 



|u| | Vu| | Au| dO < C|u| L 6|Vu| i 3|Au| i 2 

ii i 
< C|u| L 6|Vu|2 2 |V9 2 u|2 2 |Au| L 2 + C|u| L 6|Vu|2 2 |Au| 

We appeal to the following inequality for functions in M 3 (the scaling 
is 2D due to the vertical integration. For a proof, see [6]), 



o 



H 



vu| dz \f\\ g \ an < c\f\ L MlMUg\h\g\l 



1 1 

1 ! 



to obtain (by setting / = Au, g = d z u) 

j (^j |Vu|d^ |9,u||Au|dfi < C7|Vu|J 2 |9 2 u|| 2 |V9 2 u|J 2 |Au|f 2 . 

The last two integrals in the equation above are handled by (2.44) with 
k = 2. After these intermediate steps we derive from (2.76) that 

d 4 
— |Vu|^ 2 + fi\Au\ 2 L2 + u\Vd z u\ 2 L2 < -\F\ 2 L2 + 

(Jiv jJt 

C (|u|te + |^u|| 2 |Va z u|| 2 ) |Vu|i 2 + C 2 \u\ 2 L2 . (2.77) 

We appeal to the Gronwall inequality, and, with use of (2.46), (2.68), 
and (2.75), we obtain 

\Vu(-,t)\ 2 L2 + n [ |Au|| 2 ds + r/ f \Vd z u\hds < J v (t), (2.78) 
Jo Jo 

with 

J v (t) = e c ^ +J ^) t+j2 ^ [\u Q \ 2 m + C 2 U + i jf |F(-, s)\ 2 L , ds^j . 

We now prove Theorem 2.1. 

Proof. The short time existence and uniqueness of the strong solutions 
of (2.1)-(2.8) can be established as for the viscous primitive equations 
(see [14, 17]). Let u be such a strong solution corresponding to the 
initial data u with the maximal interval of existence [0, T*). If T* > T, 
then there is nothing to prove here. Let us suppose that T* < T, and 
in particular, T* < oo. Then it is clear that 

limsup ||u||#i = oo. (2.79) 

Otherwise the solution can be extended beyond T*. However the es- 
timates (2.75) and (2.78) indicate that ||u(-,i)||#i < oo for all t <T, 
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which contradicts (2.79). Hence the solution must exist for the whole 
period of [0, T). 

It remains to show the continuous dependence of the solution on the 
data, of which the uniqueness of the solution is a consequence. Let us 
assume that u 1 and u 2 be two solutions corresponding to the two sets 
of initial data and u 2 ,, respectively. We let 



u = u — u 

Then u satisfies the following equation 

~dt 



u 1 ■ V)u+ (V ■ J\ l + (u ■ V)u 2 + (V ■ |"u 



du 2 

dz 



d 2 u d 2 
+ /kxu - fiAu - v— = -fiAP Ml NU - u—P MjN u. (2.80) 



We multiply (2.80) by u and integrate by parts over Q to obtain 

- — I 
2 dt' 



±iu|£ a + ((u 1 .V)u,u) 



(V ■ j" u 1 dO^, uj + ((u • V)u 2 , u) + ((V ■ jf udO^, u) + 



Vu| 



V 



du 



dz 



fi\VP M ,Nu\ +U 



-^-Pm,nu 

dz 



We verify that 



((u 1 ■ V)u, u) + ((V • jf u 1 dO^, u) = 0, 
|((u- V)u 2 ,u)| < |Vu 2 | i2 |u|| 2 |Vu| 



3 

I 2 

lL2) 



and that 



du 2 



Dz 



(V • / ud^)— , u ) < C|u|| 2 |Vu|| 2 |^u : , L2 



\Vd 7 u 2 



L 2 ' 



(2.81) 

(2.82) 
(2.83) 

(2.84) 



Again, the last two integrals on the right-hand side of (2.81) are han- 
dled by (2.44) with k — 1. Thus we derive from (2.81) that 



dt 



u 



2 

L 2 



/i|Vu| 2 + v 



du 



dz 



< 



C (|Vu 2 | 4 L2 + |^u 2 | 2 2 |V<9 z u 2 | 2 2 ) |u| 2 2 + dlull,. (2.85) 

Thanks to the a priori estimates (2.75) and (2.78) and the Gronwall 
inequality, we have 

|u(-,t)|| 2 < e c ^(*)*+^(*))+c 1 | u (. )0 )| 2 2 . (2.86) 



22 



QINGSHAN CHEN, MAX GUNZBURGER AND XIAOMING WANG 



This shows that the solution depends on the initial data continuously. 
When u(-, 0) = u x (-, 0) - u 2 (-, 0) = 0, 

u(-,t) = 0, for all* >0. (2.87) 

This shows the uniqueness of the solution. □ 

3. A SPECTRAL- VISCOSITY MODEL FOR GEOPHYSICAL TURBULENCE 

In this section we study a model that has applications in the simu- 
lation of geophysical turbulent flows. 

3.1. The model. The 3D primitive equations with linear spectral 
eddy viscosity read 

dn , „. du 1 „ , d 2 u , „ 

— + (u ■ V)u + w— + fk x u + — Vp - /iAu - v—r - Lu = F, 

ut OZ pQ oz 

(3.1) 

dp 

^- = -Po9, 
oz 

(3.2) 

„ dw 
V-u + — = 0, 

oz 

(3.3) 



Lu = u s A(I - P m ,n)u + v&— (I - Pm,n)u. 



with 

' dz 

The other notations being the same as those in Section 2, the newly 
introduced ones n$ and v$ are the eddy closure parameters in the hor- 
izontal and vertical directions, respectively. The subscript 5 indicates 
that the parameters depend on the grid resolution. 

As for the model with partial viscosity in Section 2, we consider a 
rectangular domain f2 = M x (— H, 0), with M = (0, L x ) x (0, L y ), 
and we also consider periodic boundary conditions in both the x and y 
directions, and free-slip, non-penetration boundary conditions in the 
vertical directions. More precisely, 

u(x,y, z,t) = u(x + L x ,y,z,t), (3.4) 
Ji(x,y,z,t)=u(x,y + L v ,z,t), (3.5) 

— \ Z=0 = —\ Z= _ H = 0, (3.6) 
on On 

w\z=o = w\z=-h = 0, (3.7) 
p\z=o = Po- (3.8) 
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As a consequence of the boundary conditions taken here, the spectral 
low-pass filter Pm,n can be, and is defined as in (2.9). 

For the model (3.1)-(3.8) we obtain two results. The first is the 
global well-posedness of the model, and the second is concerned with 
the convergence of the solutions of (3.1)-(3.8) to those of the viscous 
primitive equations as fj,s and us tend to zero. The first result shall come 
as no surprise, and therefore its proof will only be briefly sketched. The 
second result will be discussed in more details. 

As usual, p and w can be expressed in terms of u. Specifically, 
integrating (3.2) from z to 0, and using the boundary condition (3.8) 
we obtain 

p(?,y,z,t)=p (x,y,t)-p gz, (3.9) 

Integrating (3.3) from z to 0, and using the boundary condition (3.7), 
we obtain 

/•O pO 

w(x,y,z,t)= / V-ud£ = V- / ud£. (3.10) 

J z J z 

Setting z = —H in (3.10), and using (3.7) again, we find 

V- / udz = 0. (3.11) 

J-H 

We substitute (3.9) and (3.10) into (3.1), and obtain a single closed 
equation for u: 

— + (u • V u + (V • / udC — + A x u + -Vp - 

ot J z dz p 



/iAu - u— - Lu = F. (3.12) 



Po 

'dz* 

The prognostic variable u satisfies the equation (3.1), the constraint 
(3.11), and the boundary conditions (3.4)-(3.6). To complete the sys- 
tem, we also require u to satisfy the following initial condition: 

u(x, y, z, 0) = u (x, y, z). (3.13) 

The diagnostic variables w and p are given by (3.10) and (3.9) respec- 
tively. 

We now specify the barotropic and baroclinic modes of (3.12). We 

let 

1 f° 

u(x,y,t) = —J u(x,y,z,t)dz, (3.14) 

which denotes the barotropic mode of the primitive variables. We also 
denote by 

u'(x, y, z, t) = u(x, y, z, t) - u(x, y, t) (3.15) 
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the baroclinic mode of the primitive variables. 

As in Section 2, we derive the equations and conditions that the 
barotropic mode u and the baroclinic mode u' must satisfy. They are 
summarized as follows. The equation and boundary conditions for u 
are 



+ (u • V)u + [(u' ■ V)u' + (V ■ u')u'] l/kxu| — Vp - A*Au = F, 

at Po 

(3.16) 
(3.17) 
(3.18) 
(3.19) 



u(x,y,t) = u(x + L x ,y,t), 
u(x,y,t) = u(x,y + L y ,t). 
V • u = 0. 



The equation, boundary conditions, and constraint for u' are 



du' 

~dt 



f° du' d 2 u' 
+ (u' ■ V)u' + (V ■ J u'dO^ l/kxu'- ^Au' - v-^ - Lu' 



Oz 



+ 



V ■ V)u + (u • V)u - ^(u' • V)u' + (V • u')u') 



F', 



W(x,y,t) = u'(x + L x ,y,t), 
u'(x,y,t) = v!(x,y + L v ,t), 



dz 
o 



du' 

dz 



z=0 

u' dz = 



z=-H 



(3.20) 
(3.21) 
(3.22) 

(3.23) 
(3.24) 



H 



3.2. Existence and uniqueness of strong solutions. The func- 
tional settings are the same as specified in Section 2. A priori estimates 
are the essential ingredients in the proof of global well-posedness of 
(3.1)-(3.8). We shall first derive the a priori estimates for u, u and u'. 

L 2 estimates 

We take inner product of (3.12) with u, integrate by parts over Q, 
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using the boundary conditions (3.4)-(3.8), we obtain 



;^, u ) + ((u- v)u,u) + ((y-J ud 0§p u 

(fk x u, u) + — (Vp , u) + yu|Vu| 2 + U 
Po 



du 



Mai V(J - P A /,at)u| 2 + u s 



d_ 

dz 



(I - P m ,jv)u 



dz 

2 



(F,u). (3.25) 



Following similar steps in Section 2, we obtain that 



|u(-,t)| 2 +/i f \Vu\ 2 ds+u [ 
Jo Jo 



-u 



"6 



d_ 

dz 
d 

— (I-P mn )m 



ds+ns / |V(/-P A /,Jv)u| 2 ds+ 



ds<ifi(f), (3.26) 



where 



K x (t) = \u \ 2 + C [ \F\ 2 ds. 

Jo 



(3.27) 



We note that Ki(t) is a non-decreasing positive function in t. This 
notion will be useful later in the paper. 

L e estimate on u' 
We rewrite (3.20) as 



— + (u'-V)u'+f V • J u'd£j ^+fkxu'-(^ 5 )Au'-(u+us] 



d 2 u' 

lh 2 



+ 



[u ■ V)u + (u • V)u' - f(u' • V)u' + (V • u')u' 

d 2 

1 dz 



¥'- 



[i s AP MiN u - u s —P MiN u. (3.28) 



Now (3.28) is in the same form as equation (2.30), and therefore the 
same techniques from Section (2) can be applied to yield the L 6 esti- 
mates on u'. We omit the details and state the results as follows. 



|u'(-,t)||e <K e 



with 



U oIl6(Q) 



c 



\L 2 (n) 



(3.29) 



ds . 
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We also have 

(li + (i8)J* J (|Vu'| 2 |uT + |V|u'| 2 | 2 |u'| 2 ) dttds + 



+ vs) 



JO, 



du 1 



dz 



'|4 



U 



— lu'l 2 

dz 



u'| 2 I dfids < K 6 (t), (3.30) 



with 



K 6 (t) = |u || e + Kl{t) / \F(-,s)\ 2 L2 ds + 



K%(t) (CKl(t) + CK x {t)t + CK^t) + d(t)) . 



Estimate |Vu|^2( M ) 

We multiply (3.16) by — Au and integrate by parts over M to obtain 



2 ~dl 



Vul^^+^IAul^^)- / /kxu-AudM 

Jm 



F-AudM / Vpo-AudM- / (u-V)u-AudM- 

m Po Jm Jm 



[ (u' • V)u' + (V • u')u' • AudM. (3.31) 
Jm 



We note that 



and 



fk x u- AudM = 0, 



M 



Vp ■ AudM = 0. 



M 



Following similar steps in the handling of the 2D Navier Stokes equa- 
tions, we obtain 



M 



(u- V)u ■ AudM < C|u|J 2(M) |VuU2 (M) |Au|J 2 



(My 



Applying the Cauchy-Schwarz and Holder inequalities one have 



M 



[u' ■ V)u' + (V ■ u')u' ■ AudM < 



C|Vu'H 2(n) 



u'| 4 lVu'| 2 dfi ) |Au 



\L 2 (M)- 
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Utilizing these estimates, and using Youngs inequality again, we derive 
from (3.31) that 



— |Vu| L 2 (M) + -|Au| i2(A/ ) < -|F| L 2 (M ) + 

C|u|l 2(M) |Vu|l 2(M) + qVu'|i 2(fi) + C f |u'| 4 |Vu'| 2 dfi. (3.32) 

Jn 

Applying the Gronwall inequality to (3.32), and using the previous 
estimates (3.26) and (3.30), we obtain 

|Vu(-,*)|| 2(M) + |y |Au|| 2(M) ds < K 2 (t), (3.33) 

with 

K 2 (t) = e CK *® (\u Q \ 2 m + C jf |F|i 2(M) ds + CK^t) + CK 6 (t)^ . 



To estimate \d z u\ L 2^ 

We multiply (3.12) by d 2 u/dz 2 and integrate by parts over Q, and 
utilizing the boundary conditions (3.4)-(3.6), we have 



d 2 u 



dz 2 



+ 



L 2 (Q) 



H & \Vd z (I - P MtN )u\ 2 L 2 {n) + u s 



d 2 (I-P M , N )u 



dz 2 
o 



/ n ((u.V)u + (v./ud Z )| 



L 2 (0) 

du\ d 2 u 



J dz 2 



' dz 2 
da (3.34) 



By Holder's inequality, 



d 2 u 

Ik 2 



< C\F\ 2 L2 



d 2 u 



L 2 
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With regard to the last integral in (3.34) we find 

du\ d 2 u 



/ n ((u.V)u + (v./ ud,) 
u- V)u4 



d_ 

dz 



udz 



dz J dz 2 
du\ du 



dz J dz 
du 



dQ 



d 2 u\ du 



j a ■ V)u + (» ■ V)ftu - (V • u)| + V • [ ud ? ^) gz 



dfl 



(d z u • V)u - (V ■ u) 



du\ du 



dz J dz 



dfl 



I -(V ■ d z u)u ■ d z u - (d z u ■ V)d z u ■ u + 2(u ■ V)d z u ■ d z udQ. 
Jn 



Hence 



u.V)u + (V.,'°ud;)|).!> 



< C / |u||V<9,u||«9,u|dQ 

Jn 

< C\u\ L e {n) \d z u\ L 3 {n )\Vd z u\ L 2 {n) 

i i 

< C\u\ L e {Q) \d z u\l 2in) \d z u\ 2 Hl(n) \Vd z u\ L 2 in) 

ii i 

< C|u|£6(n)|0*u|i 2m \u zz \ 2 L2 \Vd z u\ L 2 {Q) + C\u\ L 6 {n) \d z u\ 2 L2(n) \Vd z u 



3 



<L7|u| 4 L6 |^u|i 2 + i |u^|i 2 + ^|V^u|i 2 . 
After these intermediate steps, we derive from (3.34) that 



\dzu\ 2 L 2 {n) +Li\Vd z u\ 2 L2(n) +i 



dt 

We note that 



d 2 u 



dz 2 



L 2 (n) 



-2f, s \Vd z (I-P MtN )u\ 2 L2(Q) + 



d 2 (I-P MtN )u 



dz 2 



<C\F\ 2 L2 + C\u\Ud z u\h- (3.35) 



|u| L e = |u + u \ L e 

< \U\ L 6 + |U'] L 6 

< C|VU| L 2 + |u'| L6 . 

Therefore, by (3.29) and (3.33), we have 

|u| 4 L6 < C (K 2 (t) + Kl{t)) • 



(3.36) 
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Then (3.35), together with (3.36), gives 



d 



-fo\dzu\ L 2 (Q) +tx\Vd z u\ L 2 {a) +v 
d 2 (I-P M>N )u 2 



d 2 u 



dz 2 



L 2 (Q) 



-2fis\Vd z (I-P M:N )u\ 2 L2m + 



2v s 



dz 2 



L 2 (n) 



<C\F\ 2 L i+C{K 2 (t)+K 2 (t))\d z u\li- 



(3.37) 



An application of the Gronwall inequality to (3.37) readily gives 

\d z u{-,t)\ 2 L 2 {n) +/i / \Vd z u\ 2 L2(n) ds + 
Jo 

\-Q^\h(n)ds + 2fi s J \Vd z (I - P^,jv)u|| 2(n) ds + 



2v 8 



t 



d 2 (I - P m ,n)k i2 

dz 2 1l2 ^ 



ds < K z (t), (3.38) 



with 



K z (t) = e c{m )+ Kl(t))t A +C J* | Fj | 2 ds ^ 



To estimate |Vu|^2 

We multiply (3.12) by — Au and integrate by parts over f2, 
^1 Vu l! 2( Q) + I ((u • V)u + (v • J° u d£j ^ ■ (-Au) dQ + 

/ / k x u • (-Au) dQ + — I Vp ■ (-Au) dfi + /i|Au|| 2(Q) + 

HVd*u|| 2(Q) + ^| A g _ p M)JV ) u || 2(n) + i/j|V9 z (7 - PM,Jv)u|| 2(n) 

= (F, -Au) (3.39) 

We note that 

/ (/k x u) ■ Audfi = 0, 

Jai 

and 



Vp • (-Au) dtt 



Vp ■ / ( Au) dz dM 

M J-H 



-H / Vpo-AudM 

JM 
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It is easy to see that 



/ n ((u.V)„ + (v./ud ? )|).(A„)d n < 

C j Mu||Vu| + J |Vu|d^|<9 2 u|^ |Au|da 



Therefore we derive from (3.39) that 



+^|A(J-P MiiV )u|| 2(n) +z/ 5 |V^(J-P^ J v)u|| 2(n) < V|£ 2 + ^|Au|| 2 

+ cj (ju||Vu| + J \Vu\dz\d z u\^j |Au|da (3.40) 

By the Young's inequality and the Sobolev interpolation inequality 
(A. 4) for L 3 functions, we have 

|u||Vu||Au| dQ < C|u|i6|Vu| L s|Au| i2 

C < |u| L 6|Vu|| 2 |V9 2 u|| 2 |Au| L2 + c7|u| L c|Vu|| 2 |Au| 
We appeal to the following inequality, 

r-0 



|Vu| dz \f\\g\ dQ < cifiM^iui^igiiMlu 

n \J-h J 
to obtain (by setting / = Au, g = d z u) 

|Vu| dz) |9*u||Au| dtt < C7|Vu|| 2 |^u|| 2 |V9 z u|| 2 |Au|| 2 . 
n \J-h J 

After these steps, we infer from (3.40) that 

^|Vu|| 2 +/i|Au|i 2 + 
v\Vd z u\h + 2fJL S \A(I - P m ,n)u\ 2 l2 + 2u s \Vd z {I - P M , N )u\ 2 L2 < 

-\F\ 2 L , + C (|u|£ 8 + |£U|| 2 |VcU|y |Vu|| 2 . (3.41) 



PARTIAL AND SPECTRAL- VISCOSITY MODELS FOR GEOPHYSICAL FLOWS1 



We are ready to apply the Gronwall inequality, with use of (3.36), 
(3.38), to obtain 

|Vu(-,t)|| 2 + ^ f \Au\ 2 L2 ds + v [ |V<9 z u| 2 2 ds + 
Jo Jo 

2(j, s [ | A(J - P m ,n)u\% ds + 2v s [ | Vd t (I - Pm,jv)u|| 2 ds 
Jo Jo 

<K v (t), (3.42) 

with 

Ky{t) = e C(Kl + K!(t))t + KUt) ( | Uo |^ + 1 f | F (., s) |2 2ds 

V ^ Jo 

Now that the key estimates are in place, with an argument similar 
to that for Theorem 2.1 in Section 2, we can show 

Theorem 3.1. For a given T > 0, let F e L 2 (0, T; L 2 (fi)) ; u G V. 
Then there exists a unique strong solution of the system (3.1) -(3.8) 
which depends continuously on the initial data. 

The proof is omitted. 



3.3. Convergence of the solutions. In this section we take fis — > 0, 
ug — > in (3.12), and study the convergence of the solution of the 
system. We first rewrite (3.12) as follows: 

du 5 ,„ f° x ,^du 5 . x 1 



u 5 -V)u 5 + (V- / u a d£)— + /kx u 5 + — Wp 

Jz oz po 



dt 

- fiAu 5 - u^- - v 6 A(I - P M , N )u s -u s ^(I- P m ,n)u 5 = F. 



(3.43) 



The superscript 5 emphasizes the fact that the solution u 5 depends 
on the spectral viscosity parameters fis and u$, which themselves are 
determined by the grid resolution. The proper relation between the 
parameters fis and u$ and the grid resolution 5 is the subject of a 
separate endeavor, and will be presented elsewhere. In this work we 
focus on the behavior of the solution u,5 of system (3.43) as fis and u$ 
tend to (corresponding to the scenario when the grid becomes finer 
and finer). 
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Let u be the solution of the primitive equations without artificial 
viscosities, that is, u is the solution of 



du f° du 

— + (u- V)u+(V-y ud£)— + /kxu + 



1 <9 2 u 
— Vp - uAu - v— - = F. (3.44) 
Po oz 2 

It can be shown that the system (3.44) plus (3.4)-(3.8) has unique 
global strong solutions, provided that the initial data and the forcing 
are sufficiently smooth. See [7, 20, 21] (These works use boundary 
conditions different from ours, but the case with periodic boundary 
conditions can be handled as well.) 

To study the convergence of the solution of system (3.43), we sub- 
tract (3.44) from (3.43), let v 5 = u 5 — u, and we have 



+ (u 5 -V)v 5 + (V-y u 5 dO^ + (v 5 -V)u+(V-y v 5 d0^j + 



dt 

fkxv 5 -fiAv 5 -v^-fi s A(I-P M)N )u 5 -v s — (I-P M , N )u s = 0. 

(3.45) 

We take the inner product of (3.45) with v 5 , and integrate by parts 
over Q to obtain 

1 d 

2 dt 

J (V- v)u + (v-^Vdo|^) Vdft + 



Wr S \ 2 4. 




I V \L 2 + 









/k x v 5 ■ v 5 an + fi\ Vv* 5 ! 2 + v 



dv 5 2 



dz 



-fi s {V{I - PM, N )u d , Vv d ) - u s ( J^(J - Pm,jv)u, jU a ) . (3.46) 



We note that 



/ /kxv*-v J dQ = 0. (3.47) 
Jn 

We can also verify by integration by parts, and using the boundary 
conditions (3.4)-(3.6) for u, u 5 (and therefore for v 5 ), that 

J + (V ■ J\ 5 d^^pj -v a dft = 0. (3.48) 
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By Holder's inequality and the interpolation inequalities (A. 4) and 
(A. 5) in R 3 , we find that 



/ ((v 5 -V)u) 
Jn 


■ v 5 dfi 


< C7|Vu| L2 |v <5 | L 3| 


[ ((v 5 -V)u) 
Jn 


■ v 5 dtt 


< C|VuU 2 |v <5 || 2 | 



\Ifi, 



(3.49) 



Using Holder's inequality again we find that 




We first note that 



i\l 





\Vv s \ dz 

■H 







C|VvV(n 



|Vv 5 | 2 dz dx dy 



M J -H 
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By the integral Minkowski inequality (A. 6) and the interpolation in- 
equality (A. 2) in R 2 , we find that 



M 



du 



dz 



dz) dxdy] <C 

< C 

< C 

< c 



H \J M 



du 



dz 



dxdy 1 dz 




dz 



L 4 (M) 



L 2 (M) 



.du(-, z) 



dz 



L 2 (M) 



du 






dz 






L 2 (n) 





L 2 (n) 



Similarly, we have 



M \J-H 



v 5 | 2 dz) dxdy] <C|v d | L2 |Vv d 



L 2 ■ 



Therefore, 



du 



< c 


du 


i 




2 










dz 


L 2 (n) 


dz 


L 2 {fl) 



V • / v d d£ ) ■ v 5 dQ 



With (3.47)-(3.50), we derive from (3.46) that 

2 



r S\ 2 I Y7-r5 I 2 



\l 2 



Vv 

(3.50) 



L 2 ' 



^|v'li, + ,|VvY 



dv 5 



dz 



<(JL S \V(I-P m ,n)v*\ |Vv*| + 



d_ 

dz 



(I - P m ,n)u 









du 


i 

2 






dz 






dz 


L 2 (Q) 



,du 

dz 



r<5 I 2 I V7t r t> I 2 



L 2 {n) 



\L 2 



Vv fl |£ 2 + 



L7|Vu| L 2|v 5 || 2 |Vv 5 || 2 . (3.51 



By Young's inequality, we have 



jjlv'li. + HVv 



<5 1 2 



<2^\V(I-P M , N )u s \ 2 + 
u 1 1 



U 



d 



dz 



(I - P m ,n)u 





du 


2 










dz 


L 2 (n) 


dz 



+ |Vu|^ 2 |v 



S 1 2 



\L 2 - 

(3.52) 

We notice that the a priori estimates obtained in the previous section 
are independent of fi$ and us. We apply the Gronwall inequality to 
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(3.52), and, utilizing the a priori estimates (3.26), (3.38) and (3.42), 
we obtain 



v a (-,t)| 2 2 + / u\Vv 5 (;s)\ 2 ds + v 



t)\h + [ /i|Vv 

JO 



dv 5 , , 2 

s. 



dz 



ds 



< e m) + KUt)(^± + VA\ Kli t). (3.53) 



We have just proved 

Theorem 3.2. Let T > be given, and the other assumptions be the 
same as in Theorem 3.1. Then, as fig&LUs — > 0, 

u 8 - u ~ y/fj, s + n — ^0 in L°°(0,T;L 2 (ft)), (3.54) 

u s - u ~ y/iis + v& — ^0 in L 2 (Q,T;H\Q)). (3.55) 

Appendix A. Some inequalities 

We list here some functional inequalities that are frequently used in 
this paper. 

LP interpolation inequality 

Let ft C M 3 , and 1 < p x < p < p 2 , Pi ^ p 2 - Let u e L Pl (ft) n L P2 (ft). 
Then u G L P (Q), and 

\\u\\lp < \\u\\lpi\\u\\lp2, (A.l) 

with 

_ Pi P2 - P P2 P - Pi 

s l — j ^2 



P P2~Pl P P2 - Pi 

Ladyzhenskaya/Sobolev inequalities in M 2 

Let M C M 2 be a bounded domain with piecewise smooth boundaries. 
For each <fi e H 1 (M), the following inequalities hold: 

H0ll^(n)<C||0||I 2(n) ||0||| 1(a) , (A.2) 

(A.3) 

Ladyzhenskaya/Sobolev inequalities in 1R 3 

Let ft C M 3 be a bounded domain with piecewise smooth boundaries. 
For each e if 1 (ft), the following inequalities hold: 

»m<c\ H li^l |0| || 1(n) , (A.4) 
lhq) < C\\(j)\\ HHn) . (A.5) 
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Minkowski integral inequality (for p > 1 ) 

Let Q\ C K mi and f2 2 C M m2 be two measurable sets, with nil and 
m 2 being positive integers. Let /(£,?7) be a measurable function over 
Q 1 xQ 2 . Then 

([ (f \m,v)\dv)" at)* < [ (f \m,v)\ p <%Y d7] - (A - 6) 

The inequality (A.l) can be verified by the Holder's inequality. For 
(A.2)-(A.5) we refer to such classical texts as [1, 33]. A proof of (A. 6) 
can be found in [16]. 
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